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Abstract 

We develop a formalism for performing real space renormalization 
group transformations of the "decimation type" using low tempera- 
ture perturbation theory. This type of transformations beyond d = 1 
is highly nontrivial even for free theories. We construct such a solu- 
tion in arbitrary dimensions and develop a weak coupling perturbation 
theory for it. The method utilizes Schur formula to convert summa- 
tion over decorated lattice into summation over either original lattice 
or sublattice. We check the formalism on solvable case of 0{N) sym- 
metric Heisenberg chain. The transformation is particularly useful to 
study models undergoing phase transition at zero temperature (vari- 
ous d = 1 and d = 2 spin models, d = 2 fermionic models, d = 3,4 
nonabelian gauge models...) for which the weak coupling perturbation 
theory is a good approximation for sufficiently small lattice spacing. 
Results for one class of such spin systems, the d=2 0(N) symmetric 
spin models (A^ > 3) for decimation with scale factor r] = 2 (when 
quarter of the points is left) are given as an example. 
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1 Introduction 



The renormalization group (RG) transformations is one of the most powerful 
and frequently used conceptual as well as practical tools in statistical physics 
and quantum field theory. While conceptually the idea of combining variables 
on neighbouring sites into complexes is very simple, in practise it almost 
always turns out to be rather complicated. 

Historically, the usefulness of RG transformations was realized after the 
d=2 Ising model on triangular lattice was very elegantly solved by Niemeijer 
and Van Leeuwen [|l| using block spinning. The nonlinear RG method they 
used was peculiar to that particular system and didn't allow generalization 
to more complicated cases. For more general systems Wilson proposed to 
use the weak coupling (low temperature) perturbation theory in momentum 
space. This was first applied to scalar 0*^ models and subsequently to spin 
systems and lattice gauge theories for thinning by factor 2 [Q. For com- 
plicated systems like these with local gauge symmetries some approximate 
methods were developed like Migdal-Kadanoff approximation, variational 
RG mean- field RG ^ or block spinning using Monte Carlo numerical 
methods p, However, unlike perturbation theory, these approximations 
are uncontrollable in a sense that it is not clear how to estimate errors. 

Exact RG transformations are generally not known (exceptions are dec- 
imations in spin chains d = 1 and mentioned above very special cases in 
d = 2). Moreover after one RG transformation the resulting action contains 
generically infinite number of interaction terms, and therefore one is forced 
to make an additional approximation dropping some of them (hopefully the 
less relevant ones). RG transformations are especially useful when applied 
repeatedly. This requires self similarity of the approximate effective action 
and is justified only around fixed points. 

Note however that accurate thinning of the lattice even just by factor 
7] = 2 can greatly facilitate the study of a model by means of subsequent MC 
simulation. 

There are several types of the RG transformations. The conceptually 
simplest one is the decimation or thinning of degrees of freedom in the con- 
figuration space. Some degrees of freedom located, for example, on sites with 
at least one odd coordinate are simply integrated out. 
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Fig. 1. Decimation: full circles belong to sublattice £*, 
while empty circles x eV = C — C* denote the integrated out sites. 



Example is given on Fig. 1 on which spins at empty circle points are 
integrated out. 

<f>{X):X&C* <t,{xy.xeC-'L 4>(X):X&C* 

Here and in what follows points of the coarse lattice are denoted by capital 
letters. The resulting effective action A'^'^'^ contains generally interactions 
of any range. Here C is the original d dimensional lattice while C* is a 
sublattice. Note that remaining variables are all the old variables. This is 
not the case for the so called block spin transformations. One defines a linear 
or a nonlinear combination of the variables on £, the block spin: 

0(X) = /[X,0(x)] (2) 

For example, for the 0{N) classical spins one can define |^ 

s'^(x)= y: s''{x)/\ e ^"wi (3) 

blockX blockX 
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This combination is highly ambiguous and success of the transformation crit- 
ically depends on it's choice. The main problem is that it is extremely dif- 
ficult in practise to perform such a transformation even perturbatively. The 
relations like eq. are very nonlinear and even singular [|1^. 

Another type of RG transformations, used especially extensively in field 



theory, is the momentum space RG g, |ri|] . One defines the momentum space 
variables 

(f){p) = 1/{27tY J rf'^xe*f"0(x) (4) 

Now one performs integration over high frequence modes (strictly speak- 
ing chopping the Brillouin zone, but more often the approximate spherically 



symmetric momentum cutoff A is utilized [y, [12|). This type of RG trans- 
formations, while convenient for the 0^ model, turns out to be especially 
inconvenient for constrained systems like the 0{N) symmetric spin models. 
The reason is following. While generally in x - space the constraint are local, 
for example 

S''{x)S''{x) = 1 (5) 

in p - space it becomes a convolution. What does it mean now high frequency 
physical modes? The constraint mixes between low and high frequencies. 
Since most systems of interest belong to this class one has to circumvent the 
difficulty. One way is to solve the constraint and make the momentum space 
RG for physical quantities only. Then the mode integrated effective action 
contains generally " non-covariant terms". The original global symmetry is 
lost since the high frequency modes do not constitute an 0{N) symmetric 
set. Problems are more acute with local gauge symmetries. In practice this 
type of thinning out of degrees of freedom is often used for the demonstration 
purposes only and very rarely the actual calculations. 

Decimation are extremely difficult to perform even in free theory in more 
then one dimension (see for example ll3|). It might sound surprising that 



something is difficult in free theory since all the integrals are Gaussian and 
"doable in principle". Of course it is still a Gaussian integral, but a very 
complicated one. Let us consider a free massless boson nearest neighbours 
action 
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where the lattice Laplacian is defined by 



If one tries to integrate out a point 0(0, ...,0) the Gaussian integral in- 
volves all its 2d nearest neighbours. 



/ dm 



exp 



2 V 



(2^0(0)^ - 20(0) $:(</.(/.) + </.(-/.) 



= exp 



2d 



(7) 



This is very simple. However when trying subsequently to integrate another 
point, say (2,0...,0), all the previous point's neighbours enter the Gaussian 
integral and so on. The Gaussian integration requires inverting increasingly 
larger matrices. Since we have to integrate out all the points not belonging to 
the sublattice, some other methods are required. An exception is the d = 1 
case. Here the size of the matrix does not grow: integration of a point leads 
just to interactions of the neighbouring remaining points. This is the reason 
why it is possible in many cases to explicitly find decimations va. d— 1. 

In this paper we perform the decimation for multidimensional free the- 
ories. The result does not coincides with the naive continuum limit even 
in the limit of large 77 = A/ a. This is discussed in section 2. Then using 
this result we develop in section 3 general perturbative formalism for weakly 
interacting models. It is apphed in section 4, 5, 6 to the 0{N) symmetric 
nearest neighbours interaction spin model (the nonlinear cr- model). In sec- 
tion 4 wc derive the general diagrammatic technique for such a models. In 
section 5 the solvable d = 1 model is considered and results compared with 
the usual perturbative ones, while in 6 the two dimensional asymptotically 
free model is studied. The perturbative method generally is better suited for 
asymptotically free (=phase transition at T" = 0) models like this one. This 
is because on fine lattices coupling becomes small. Then the RG transformed 
model can be investigated, say by the MC method, on coarser lattice. Even 
simplest decimation with rj = 2 reducing the number of points by factor 4 
greatly simplifies the numerical work. 
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A coarse grained effective action is generally obtained as a series in " close- 
ness" of the interacting spins: the nearest neighbours, next to nearest etc. 
In order to make any practical calculations possible, one has to truncate it 
at some point. We restrict our consideration here to the fourth order terms 
in fields and up to the fourth derivatives. 

We conclude in section 7 by discussing complexity of such calculations 
and some of their uses. 



2 Decimation of free fields 

Let us start with free boson theory on the lattice. Effective action, after 
decimation with parameter rj generally has a form: 

A^m)] = ^ E 'PiX)MX - Y)ct>{Y) (8) 

where bold letters denote sublattice functions. Of course it is quadratic in 
Let us now perform Fourier transforms of the original and sublattice fields 

'^W=(^i:e'^^W- (9) 
This convention fixes the Fourier transforms of propagators: 

G{x) =< <P{x)m ^ £^\'ke^'^^G{k), 

G{k) = a'^Y.^~'^'"'G^^)^ (10) 

G(ir)=Ee-^^^GW, (11) 

X 
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and inverse propagators 

m = ^ = ^Ee-''^"A(a;) = A^sin2(aV2), (12) 
A(ii') = ^e^^^A(X) (13) 

Due to translation invariance, we can invert these in momentum space to 
obtain propagators on the lattice and sublattice correspondingly: 



1 f-K/a 
^ ' (27rYJ-n/a 



-ikxa 



(2x)<'i-,/a 4E>iin'jV2) 

The two models should result in equivalent correlator between two sub- 
lattice points: and X: G{riX) = G{X). This leads to the following relation 
between the Fourier transforms: 

G(K) = 1/A(K) = -i-y e-^^^ T^" dV ^ , , , (15) 

^ ^ ' ^ ^ {27TYjt J-/- E;.4sin2(aV2) ^ ^ 

Summation over X results in sum over S functions 

^^^) = 7^/ . ^ ^((fc-K + 27rn),) (16) 
{2nyj-^/a „fr'oo E^4sm(aV2) 

which are used to perform the momentum integrations: 

^^^^ ^ „5i E,4sin^(a(i^ + 27rn),/2) ^^^^ 

Limits of summation in the last expression follow from the different sizes of 
Brillouin zone for two lattices (see Fig. 2). Note that due to periodicity the 
limits of summation in n can be shifted by the period rj. 
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In (i = 1 we recover the previous result since the sum is doable |T^. In 
d > 1, the summation over one of the variables, ni can be performed similarly, 
but the remaining summations should be done numerically. In particular for 




Fig. 2. Brillouin zones for original lattice (k) and sublattice (K). Lines 
correspond to k = K + 27m. 

d = 2 we have the propagator: 

= 1 V ^i^h(^c.)csch(c.) ^_ 

^ nf^i 2 (-1 + cosh(r7a) + 2sin2(f )) 



where 



2 K2 + 27^2, 



a = arccoshfl + 2sin ( )) 

27] 



For large rj the Euchdean invariance is restored, G{Ki, K2) = G{y + Kl) 
and numerical calculations show it can be fitted by 

G{K) = + —\og{r]) + 0M876 + 0.003022 + 0[{Ky] (19) 

K'^ 277 

(see Fig. 3) with an accuracy better 1 percent in all the Brillouin zone. The 
first term is the continuum propagator. Note that the decimated propagator 
even for large rj does not coincides with the naive continuum limit. 

The contact constant term with logarithmic dependence of rj is typical for 
d = 2 and is nothing else but the bubble integral. The polynomial coefficients 



8 



are very small and almost coincide with the Loran 




expansion of the propagator around K = 0. For finite rj the symmetry 
remains of course just the discrete subgroup of the rotations. 

In higher dimensions similar expressions can be written. Similar proce- 
dure can be extended to free fermion fields (see Appendix A). Especially 
interesting aspect of this is the species doubling |T^ . 

In the simplest case of one dimensional massless boson field we can ex- 
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plicitly integrate out all the odd points since the integrals do not intertwine: 



n 



b2x+i exp 



exp 



4>2x+l ~ i.'t>2x + 4>2x+2)<t>2x+l 

X 



2 (2a) 



>xvx+i 



(20) 



We therefore obtain the original form with twice lattice spacing: j32 = /5/2, 
where /3 is an inverse temperature. The action is a perfect one [|l6l. For 
arbitrary t] = a/A we get similar results, 

f3, = P/2\ 

In higher dimensions we still can perform the decimation using momen- 
tum space in the intermediate steps to diagonalize the matrices. 



3 General diagrammatic method for evalua- 
tion of the decimated action 

Now we would like to build a perturbation theory for the decimation-type 
RG transformations in the interacting case. For concreteness we discuss the 
lattice 0^ model 

X 

where Vcpx = {4>x — 4>x-i)- Low temperature (weak coupling) perturbation 
theory for this model can be represented via Feynman diagrams including 
propagator and the four vertex. In momentum space RG, when the high fre- 
quency modes from A to A' = A/rj are integrated out, the resulting effective 
action on the scale A' has a general form 

OO 1 

^''['^J = E TTTIT / r(2'^)(a:i, .., X2.)0.,...0.,„ (22) 

[2ny. Jxi,...X2n 

The coefficient functions F*^^"-* are sums of all the one particle irreducible 
contributions with 2n ends. The external momenta are all below A' while 



:(V0.)' + 



m 



4! 



(21) 
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all the integrated internal momenta are between A' and A [0. Since we 
will significantly modify the procedure in x space let us briefly outline the p 
space diagrammatics for RG. This is most easily done if original vertices and 
propagators are split into several pieces. The vertex decomposes into: the 
vertex connecting just high momenta modes (Fig. 4(c)), only low momenta 
modes (Fig. 4(d)) and mixing the two (Fig. 4(e), 4(f), 4(g)). 



b 



e t g 

Fig. 4. Momentum space RG propagators (a, b) and vertices (c, d, e, f, 
for (j)'^ model. Low momenta are indicated by bold lines. 



Propagators are decomposed analogously into low and high frequency 
parts: G{k) = e{\k\ - A')G{k) + ^(A' - \K\)G{K). Note that the coupling 
between the modes is by means of the vertex only. This will be completely 
different in real space RG. The integral over high frequencies (denoted (px) 



4'k 



exp 



k>A' 



(23) 
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with 



>K\ 



A 
4! 



•m(P{-k-l-m) 



(|fc|,|i|,H)>A' 
\ |fc+i+m|>A' 



+ 



/ 



(t>K(t>l(t>m(t>{-K-l-m) + 



|i^|<A',(K|,|m|)>A' 
\k+l+m\>A.' 



(|X|,|L|)<A',|m|>A' 
|fc+Z+m|>A' 



(t>K^L(t>m^{-K-L-m) 
\ 



+ j <t>K(t>L(t>M(t>{-k-i-m) + permutations 



(|KMLMM|)<A' 
|fc+«+m|>A' 



(24) 



apart from "classical" parts independent of (pk is exponent of the vacuum 
energy of the high frequency theory with (px playing a role of the external 
sources. This is the sum of all the vacuum diagrams in this theory. However, 
as we remarked before the lines do not connect low to high frequency modes 
and consequently all the one particle reducible diagrams vanish. 

For the real space RG the perturbation theory can be built in a similar 
way. The fields (j)x-, X & C* will be treated as "external sources", while all 
the internal points will belong toV = C — C*. For 0^ model this means the 
following decomposition (Fig. 5). Action is divided into three parts: 



-o 



-o 



a 



-o- 



Fig. 5. Real space RG (decimation) propagators (a, b, c) and vertices (d, e) 
for 0^ model. Full circles belong to sublattice (external fields), 
while empty denote "internal fields" (x e £ — £*). 
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"classical action" of the "external" field 4>x (Fig- 5(a), 5(d)), 



X 



„ Y.<t>xAxYcl>Y + Af^t[<l>xl (25) 



cross-term (Fig. 5(b)) 



0x] = -a("-') E </'xA(X, x)0, = -a('^-2) E (26) 

x,X x,X 

with "external legs" 

A(X, x) = E(5x-.+M + Sx-x-f.) (27) 

and an internal part for which all the vertices belong to "D = £ — £* (" deco- 
rated" model) (Fig. 5(c), 5(e)): 

M<Px] = E <t>M^^ y)<t>y + ^ E '^t (28) 

x,y X 

Note that unlike momentum space RG, here external fields (px are coupled 
to internal part only via derivative couplings like off-diagonal part of propa- 
gator (Fig. 5(b)) or derivative interaction in nonlinear cr-model (see the next 
Sections). All the local vertices will completely decouple into internal (Fig. 
5(c)) and external (Fig. 5(d)). 

Integration out all the fields (j)^ will lead to an effective action for the 
fields 4>x on sublattice of the form: 

A[^X]= E 7^i^('"H^l,-,^2n)0X,...0X.„, (29) 

where the coefficient functions H^'^'^\Xi^ ..^X^n)-, contrary the momentum 
space RG, are sums of all the connected contributions with 2n ends. These 
connected functions does not degenerate into one particle irreducible. 

To integrate over field (px perturbatively, we need to find its propagator 
which is matrix inverse to D. At the same time, all we know explicitly is 
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the full Laplacian A and original propagator G = A~^. Moreover, since 
V does not constitute a sublattice, it is impossible to make use of Fourier 
analysis on it. Therefore it is useful to represent all the summations over V 
via summations over C and C*. This can be done using following algebraic 
trick 0. Matrix A as well as matrix G can be decomposed into following blocs 

A=(^, ;^), (30) 

G=(:. 2). (31) 

where the infinite-dimensional matrices A, B, D, a, b, d are defined as follows. 



A, B, D defined by the quadratic part of action eqs. p5| , |27| , |28|) (Fig. 5(a), 
5(b), 5(c)). The matrices d,b are the usual propagator matrices (Fig. 5(b), 
5(c)) and a is the inverse propagator between the points of sublattice, that 
is an expression inverse to G{X). Now we can invert D using the fact that 
matrices A and G are inverse to each other. This implies a set of algebraic 
relations for their submatrices: 

Aa + Bh* = 1, Ab + Bd = 0, B^a + Db^ = 0, B% + Dd = 1 

and, after straightforward transformations we obtain an expression for the 
"internal" propagator: 

D-^ = d- b*ab (32) 
or, returning to previous notations, 

Dxy = Gxy - GxxGxyGYy (33) 
X,Y 



Here and in what follows bold line denotes free inverse decimated propa- 
gator G^^y, while thin line corresponds to propagator of the original theory 
Gxy Using this representation of D~y as an internal line, we can now begin 
to build the perturbation theory. We would like to stress out that using of 

^Somewhat similar technique based on the Schur formula was developed for quasiperi- 
odic systems on the octagonal lattice . 
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^xy ( p3D enables us to extend this summation over all the original lattice 



C rather then "decorated" subset T). Indeed, one can see that this expression 
is equal to zero when at least one of the points x, y in D~y belongs to C* (Fig. 
6). Therefore, any additional contribution due to this extension is equal to 
zero as well. 





Fig. 6. An additional contribution from extended summation. 
Small circles denote the points of the sublattice C* . 

Calculation of an n-point function H^^\Xi, in effective action eq. 

(|29|) for the coarse-grained field (px will be as follows. 

All connected diagrams with n end points Xi,..,X„, V vertices and / 
internal lines in real space are drawn with following components: 

a) all vertices are situated at the points x E C and to every vertex at 
point Xi corresponds summation J2xiJ 

b) Axx are assigned to external ends X and 

c) the internal lines D~y are represented via eq. (|55D . This representation 
splits each diagram into 2^ subdiagrams and each of these subdiagrams should 
be calculated separately. This calculation includes summation over all the 
internal points Xi on the fine grained lattice and over the internal sublattice 
points Y [end points of the inverse decimated propagator G~^, see eq. (^Sf )]. 



Finally the "classical contribution" to coefficient function should be added. 

For the sake of simplicity let us discuss calculations of the decimation 
diagrams on the concrete example. Namely, we will consider one of the 
contributions to four point function H^^^ in 0^ model (Fig. 7). 

With using of the propagator D~y, original diagram splits into five sub- 
diagrams (Fig. 7(a), 7(b), 7(c), 7(d), 7(e)). Typical subdiagram here (for 
instance, subdiagram 7(d)) can be written as 
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X! ^^1 yiGyi nGyj^y^Gyj z^X2 ViGy.^ y^Gy^YaGya z^X^ y^Gx^ z^X^ yfix^ z- 
yi,--,yi,z 

Yi,..Y4 

(34) 

As usual, in practical calculations it is very convenient to employ the Fourier 
transformed functions at the intermediate steps. In this way we will operate 
with vertices, legs 

A{k) = 2a'^-2 cos{k^a), (35) 

propagators G{k) and inverse decimated propagator G^^{K) eq. (0). How- 
ever, it is not convenient to perform Fourier transform of expression eq. (p^) 
immediately, because it contains functions defined on different lattices and, 
therefore, obeying different transformation rules (cf. eqs. (|TU|,|T2|,|TT],|T3|)). 

Instead, we can use the fact that this expression breaks into blocks where 
all internal points lie on the fine grained lattice C and sublattice points enter 
only as ends. These blocks are connected by G~^{Yi, Yj). In terms of such a 
"/^-connected" parts, diagram eq. (|34D has the form: 

= Yl ^1(^1, yi)MX2, Y2)F,{X,, X4, F3, n)Gy^V,Gy3V, (36) 
YU..Y4 

with blocks 

Fi(Xi, Y,)=J2^XAyAGy,n, (37) 

yi 

^2(^2, Y2) = Y ViGy^ Y3 (38) 

y2 

and 

Fsi^X^, X4, Y3, Y4) = Y, Gy2 zGya z^Xi ysGx-i z^X4y4,Gx4 z- (39) 

ys,y4,z 

Now we can proceed as follows. First we will Fourier transform each block 
in eq. (|36|) with respect to its internal points. This will result in diagram G^^^ 
as function on sublattice, and we can apply to it the Fourier transformation 
rules for the sublattice eq. (|Tll , p!3D . First step does not differ from the usual 
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perturbative lattice calculations, and resulting functions read, for example, 
as 

Fi(Xi, Fi) = r " e*(^i-^i)^A(A;)G(A;). (40) 

J—ir/a 




Fig. 7. Four point contribution to the 0"^ effective action 

This, however is not the case when the remaining Fourier transforms and 
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summations over internal sublattice points are performed. Due to the differ- 
ence between the Brillouin zones result will be sums rather then monomials. 
For instance, for the block Fi we will obtain: 



all nu=l 



r/2V 4sin2^:^^^i±^^ ^ ^ 



2r] I 

Such an expressions we will call "decimated" and denote by [|"expr"|]: 

[|/(X,L,...)|]= /(ir + 27rn;,,L + 27rn^,...) (42) 

Notice that the decimated function does not possess original translation in- 
variance (with the period In jd) but instead turns out to be In jijia) periodic. 
Eventually, diagram G''^ takes the following form: 

= [\F,{K)\]G'\K)[\F,{P)\]G-\P)mK, P,Q,-K-P- Q)\]. (43) 

Another simple example, diagrammatic calculation of the d = 1 free boson 
effective action, is given in Appendix B. 

Note that described algorithm is applicable not only to the tree level 
contributions. When the representation of the propagator eq (|33D is used, 
some of the loop diagrams will contain G~^. If this is the case, one should 
calculate the corresponding decimated functions, and only then perform the 
loop integration. We will encounter such diagrams in the next section. 



4 0{N) symmetric classical spins: weak cou- 
pling expansion 

In this section we apply the formalism developed in the previous section 
to the 0{N) symmetric nonlinear a model. In d dimensions, this model is 
described by action 

A = -^^5'^n^^ (44) 
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where S"' {a = 1,..,N) is 0{N) vector normalized on unity, S"^ = 1, g is 
the couphng constant or temperature and □ is the lattice Laplacian. The 
partition function of this model is given by the path integral 



Y[dS:S{Sl -l)exp 



— ys^as" 



(45) 



This is an example of a theory with constraints. To develop a perturbation 
theory for this model, it is convenient to re-express it in terms of the un- 
constrained fields |jl8| vr*, i = 1,..,N — 1. For this purpose one can solve 



constraint for , obtaining = {^cS-'^ ^^/^vr*, y'l — ga^'^^'^^n'^), and then 
the partition function in terms of "pions" vr* will have the form: 



n 



dnl 



exp 



9 X / _ 



(46) 



This last expression when used perturbatively gives rise to the infinite set 
of vertices both local (coming from the exponentiated and expanded measure) 
and via derivative couplings originating from the part of the action. For 
our present purpose it is sufficient to restrict diagrammatic to the order g"^ 
in coupling constant g. To this order, partition function takes the form: 



Z = JlidK exp 



Jd-2) 



yinW + g 



Jd-2) 



+9' 



a 



2 (d-2) 



IT j Uvr + 



Jd-2) 



(vrY +... 



(47) 



From this expression, the basic diagrams will be: massless propagator 



G{k) 



E^4sin^(afc^/2) 



and vertices (Fig. 8) 
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(a) (b) 




Fig. 8. First-order (a, b) and second-order (c,d) vertices of cj-model. The curly 
line stands for a lattice Laplacian, the broken line - for a 5 function. 



Now we can perform decimation perturbatively. First of all, one can 
see that the "classical" part of the action comes from the measure and the 

cross term ^ J2xx \J^~ 9 o^'^~^^7r|^Axx'\/l — 9 a^'^'^^n^. Indeed, other classi- 
cal terms are: 



1 

9' 



The same cancellation takes place, of course, also for the diagonal terms of 
the "internal" part of the action. At the same time the cross term expands 
as 
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9 xx 



Xx 



Ad-2) 



+ 9 



E 

X 



Q2(d-2) 



2 

Q2(d-2) 



,(d-2) 



^(cZ-2) 2 y ■ ' / 2\2 



^2{d-2) 



X 



2 \2 



+E 



(48) 



9 a 



2{d-2) 



Y^nj,Axx7rl + ... 



Xx 



X X 



\7^XMX-X) 

(a) 




(b) 





^7ri5(x - y)A{X - x) ^.7:Wy5{X - Y)A{X - x) 

(c) (d) 
Fig. 9. The lowest order cr-model sources. 



Terms belonging to the sublattice C* contribute to the classical action 
while the terms lying on T> complete the remaining off-diagonal internal part 
to the usual lattice action. 
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The basic diagrammatic elements are: free massless propagator G{x — 
y), vertices, coinciding with the usual cr-model vertices (Fig. 8), sources: 
external leg 2 a^'^''^^A{X — x) (Fig. 9(a)) and cross interaction (Fig. 9(b), 
9(c), 9(d)) and "classical" terms 



(d-2)2 3(^^f^^/2 \2 , 9 O'^ ■ J2 



(d-2) 



4 V ^/ ■ 2 



(49) 



Notice, that unlike local theories like 0^, in cr-model the "classical effec- 
tive action" contains infinite series of such decoupled terms. Performing the 
Fourier transform according the rules Sect. 2, we obtain the corresponding 
functions in momentum space: propagator 

GW = ^sin^(^), 
vertices (Fig. 8(a), 8(b), 8(c), 8(d)) 

{d-2) rp.{d-2) h. 2 „2(d-2) 



„2 _3((i-2) „ h, 

16 ^ 2 ^' 

and sources (Fig. 9(a), 9(b), 9(c), 9(d)) 



2 2{d-2) . 

^—^ 7r^7ri(5*^(5'='5™"[| A(K + L), 

^\ 'k<^%7^q5'^5^'5^^[\A{K + L + P + Q), 



where 



A(A;) = 2 a(''-^)cos(aA:), 



and "half-decimation" [|A means that once the " L-connected" part will be 
formed, it should be decimated. The classical part remains, of course, con- 
stant. 
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5 Soluble model: d = 1 Heisenberg chain 



Let us consider first the simplest nontrivial example: d — 1 a model and 

decimation with 77 = 2. We start from the two point function. The tree level 
contribution is, obviously, that of the free theory (see Fig. 17 in Appendix 
B), in this dimension 

Hi'\K) --{I- mK)g{K)\] + [\A{K)G{KWG-\K)^ . (50) 

This contribution can be calculated analytically, and one can see that it 
reduces to 

Hi'\K) = G-\K) = = 4 sin(|). (51) 

This corresponds to the fact that the free massless bosonic action in 0? = 1 
is perfect. 

One loop contribution to propagator consists of three classes of diagrams: 
diagrams coming from the measure, bubble diagrams and the self energy part. 
Measure gives the contribution shown on Fig. 10, which reduces to 

^«W = -^(l + -cos^(^))- (52) 
a ria 2 



Hg\K) = 00 + cr>_/-^ • ^-^^ 



— 2 cr>_^~^ e e • 



+ crv>~^ e e • e e ^-^^-jn^ 



Fig. 10. Contribution from path integral measure to one loop effective action. 
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Bubble diagrams are shown on Fig. 11. They give the contribution 



g{N-l) . 2(VaK, 
-sm ( 



2a 



(53) 



Notice here that although these integrals can seem to be IR divergent, in 
fact this is not the case: all the divergences cancel. The cancellation is due 
to the fact that in decimation only short distance effects are involved; long 
range "tails" remain unaffected. 



+ 





- 2 



V0 



CTV^ ' o o 




c d 

Fig. 11. One loop bubble diagrams. 



+ 



V0 CTV^ ' O O I O O 

e 




Last group, self energy diagrams, give the following contributions depicted 
on Fig. 12. This contribution reduces to 



(54) 



One can see that again all the infrared divergences canceled. 
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Crx_> 1 1 ^~Xj3 — 1 e — © ' ''~\J3 

(a) (6) 



2 ( o o ' ^ '"v© — cf~^ ' o o ' — e — © — ^ '"v© ) 

(c) id) 



+ ( cTX^ ' o o ' ' o o ^ '~VQ — cr\^ ' o o ' — e — © — ' o o ) 

(e) (/) 
Fig. 12. Self energy-type contribution to one loop effective action 

The total one loop quadratic part of an effective action with restored 
coupling constant therefore takes the following simple form: 

where □(X) is understood as lattice operator. 
Our next step will be the four-point function. 

Let us consider next terms in an expansion of the effective action. They 
contain four fields and up to four derivatives [see eq. (^7|)]. At the tree level, 
there are three contribution to these terms. First is a "classical" term A'^ 
eq. (^). Two other contributions are given by diagrams Fig. 13, 14. 

Again, using the representation eq. (PSI) splits every diagram into several. 
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(a) (b) (b) 

Fig. 13. Fourth-order contributions including internal vertex Fig. 8(b). 
M,k 



K,i K + L. 

o. 




K,iooM,k 




L,j ooM,k 
K+M ) L+M 




(a) 



(b) 



(c) 



Fig. 14. The fourth-order contributions coming from the source Fig. 9(b). 

We calculate the coefficient function H*^^) as follows: an effective action 
as decimated one is given in momentum space by diagrams (Fig. 13, 14) plus 
classical contribution: 

-AehM= [ {^'k<){^m^'-k-l-m)^h\%{K,L,M) 



K,L,M 

{'^k'^l){'^m'^-k-l-m) 



4a 



-I- (diagram Fig. 13) -|- (diagram Fig. 14) 

(56) 
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where 

(diagram Fig. 13) = -(^/4!) \([\A{K)A{L)A{M) 



X A{-K -L- M)G{K)G{L)G{M)G{-K -L- M)G-\K + L)|] (57) 

+ (terms due to (|55D)) + (permutations)] , 

(diagram Fig.l4) = (2c//4!) [([\A{K)A{L)A{K + L)G{K)G{L)\] (58) 

+ (terms due to (|33|) )) + (permutations)] . 

On the other hand, this term in an effective action on the coarse grained 
lattice in one dimension could have the only form: 

-4?/K] = - / ^H;;L(ir,L,M)vri,7ri7rW^_^ (59) 



K,L,M 

with 



H%{K, L, M) = H'{K, L,M) + c H^.^iK, L, M). (60) 



Here the first term, 
^].ki{K, L, M) = (^/4!) (%4i(-4 sm\{K + L)/2)) + (permutations) 



was determined in previous subsections (it comes from expansion of the co- 
variant action in terms of pions), and the second term is 

^l,AK, L, M) = -ig/A\) (5,/,,(4 sin2(ir/2))(4 sin2(M/2)) 

+ (permutations)) . 

To find the remaining coefficient c, we equate the coefficient functions (= 
fourth functional derivatives in fields vr)^) in expressions eqs. (56,59) at some 
definite momentum configuration, for example, for K = — L = M = P (the 
case of "back-to-back scattering") and for definite flavour indices {i = j = 
k = I = 1, for instance): 

H[tUP, -P, P) = Hffii(P, -P, P). (61) 
Then the coefficient c is determined by the linear equation 
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dP^ I I dP^ 

/ p=o ^ 



P=0 



dp^ 

For this configuration: 

-P, P)=9 (-12 - 12cos^(P/2) + 2 (5 + 6 cos(P) + cos(2 P))) , 



(63) 

HStii(P, -P, = ^ (-4 sin2(P) - 384 c sin^(P/2)) . (64) 



Equation eq. (|6^) gives then: 

c = 1/32. (65) 



Notice, that at the first sight the decimated expression, eq. (|56D, should 
behave as 1/P^ (one can check that using the representation (^) cancels 
leading term 1/P^ even without decimation). However, due to decimation 
procedure all the terms up to constant cancel, and in fact hIjIi{K, L, M) 
begins from the second order terms. Moreover, equality of the P^ terms in 
Hnn{P, —P, P) and H^^\]^(P, — P, P) can serve as some additional consis- 
tency check. 

The one loop contributions to the four derivative terms are shown in 
Figs. 18, 19 and 20 (Appendix D). For the sake of simplicity, we restrict our 
consideration to the leading 1/A^ contribution. These diagrams are calculated 
analogously to the one loop two derivatives terms, and give the following 
contribution to the effective action: 

'-^^^^Ui-r-l-'^-'). (66) 

Again, all the potential IR divergences cancel due to the decimation proce- 
dure. 

To calculate the four point function, it is more convenient to employ 
somewhat different method. Namely, instead of direct calculation of the 
effective (decimated) action, we will use a matching approach. The clue 
to this approach is that decimation does not change field variables, so that 
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their correlators between points of sublattice £*, calculated in both original 
(fine grained) and effective (coarse grained) models should coincide. We will 
therefore calculate correlators of fields in points Xi,X2, e L* starting 

from different scales, a and A in terms of the lattice spacing, and after this 
we will require matching conditions between these functions to be fulfilled. 
For practical calculations this means that wc should compare the amplitudes 
under considerations for two cases: decimated with parameter r] = Aja for 
the original lattice model and the usual lattice amplitudes calculated from 
the " phenomenological" Lagrangian including irrelevant operators with as 
yet free coefficients: 

G(")(^^i,^^2,..,r/X„) = G(")(^i,^2,..,^n). (67) 

In momentum space this condition takes the form 

[|G(-)(p,A;,..)|]|(,=p,,=^„,.) = G(")(P,X,..). (68) 

In fact, the usual approach when the lattice quantities are compared to the 
continual ones is nothing but the matching for particular value A — > oo, 
although its physical meaning is not as transparent as for the decimation 
RG. Strictly speaking, any amplitude is given as a power series in momenta, 
so that we should truncate this series at some point and consider truncated 
effective Lagrangian with finite number of irrelevant operators rather than 
exact one. Matching conditions will then fix these coefficients. Here we will 
restrict our consideration to the four derivative terms in an effective action. 
In d = 1 (7-model, there could be only one such term, [Sx^SyY. Thus, our 
"phenomenological Lagrangian" has the form: 

L = U^-{Sxn{X-Y)Syf, (69) 
9 

where Lq is the usual cr-model Lagrangian. 

Because there is only one arbitrary coefficient in a decimated Lagrangian, 

we need the four point function calculated for one particular configuration of 
external momenta, for example [|G'-^],^(P, P, — P)|] • By four point function we 
mean correlator of unconstrained fields, tt*, rather than constrained fields S"" . 
Then the second configuration (back-to-back scattering, [|G-^],;(P, — P, P)|], 
for instance) can be used as a consistency check. 
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On the tree level, this amplitude as calculated from the scale a is given 
by three contributions: 

[\G^^l,i{K,L,M)\] = S,M\G^i\K,L,M)\] + S,kSji[\G^\K,L,M)\] 

+S,i6,k[\G'i\K,L, M)\] 

To perform matching, it is sufficient to consider the term 5jj5fc«[|G^i^'*(p)|] only. 
This term is given by diagram Fig. 15 and for configuration K = L = —M = 
P has the following expansion: 



(70) 



K, io 



5ij5M[\G't\K,L,M)\] 




M, k 



Fig. 15. Fine grained four point function 



The same correlator but calculated as a sublattice quantity, has the fol- 
lowing form (Fig. 16): 



5ijSkiG^^\K, L,M) = 
K, i ^ M, k K, 



+ 2 c 





M, k 



+ permutations ) 



L, j / ^ I L, f ^ -K - L- M, I 

Fig. 16. Four point function in effective theory. Strokes correspond to the Laplacian. 
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G^^\K,L,M) = -G{K)G{L)G-\K + L)G{M)G{-K - L - M) 

-2 c [G{K)G{M) + G{K)G{-K - L - M) 

+G(L)G(M) + G(L)G(-ir-L-M)], (71) 
and expands a.t K = L = —M = P as 



G(4)(p) 



4 

P6 



p4 



+ ... 



(72) 



Comparing these two expressions, eqs. (|70|J7^) , one can find c = 1/32. Thus, 
the only tree level four- derivative term in the decimated effective action in 
d = 1 is 



X 



In fact, in the case rj = 2 the problem becomes especially simple. In this 
case the partition function has the following form: 



Yl dnl exp 



X 



X,x 



X 



^nlAxxK 

X,x 



\ " _2 



2a ^ 



TT 



X 



2a 4" ^ 



1 



4)^-EK^)^ + ^E(- 

X X X,x 



'^X) ^XxT^^ 



W-I^xxirrlf \ + 

^ X,x 



X 



(73) 



Therefore, for t] = 2 an internal line reduces to the contact term, □ 



xy 



-26. 



xy 



for (x, y) G V. This means that both an internal propagator and all the 
internal vertices are local, G{k) = a? ets (the problem becomes classical). 
The only remaining non-local terms are sources. Moreover, in this case T) is 
a lattice itself and thus there is no need to use the representation eq. (^) 
(it becomes trivial as the only non zero contribution comes from the basic 
diagram without any replacement). The loops shrink into the points and 
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therefore there is no loop integrals. Due to all these simplifications, all the 
calculations can be done immediately in the real space. 

We can compare our perturbative RG results with the perturbative effec- 
tive action obtained from the partition function eq. ([731). One can see that 
these expressions coincide. 

Effective action obtained here differs from the exact effective action for 
the d = 1 a model [^. This difference is not surprising, however. Indeed, 
an exact decimation takes into account both perturbative and instanton-like 
nonperturbative configurations while here we restrict our consideration to 
the perturbative contribution only. 



6 Two dimensional cr-model 

The real space RG discussed in previous sections is not restricted to the d = 1 
models only, but is immediately generalized to the higher dimensional theo- 
ries. In this case calculations become much more cumbersome but otherwise 
all the technique remains unchanged. Here we would like to consider one of 
the most popular two dimensional theories - d = 2 nonlinear 0{N) cx-model. 

Diagrammatics of the model under consideration coincides with one di- 
mensional case. Therefore, we can immediately use diagrams described above 
and the formalism from Section 3 to calculate an effective action. The only 
differences are disappearance of factors a^'^"^^ from the action and, of course, 
two dimensional sums in place of one dimensional. In d = 2, however, action 
is not perfect anymore and thus an exact analytical results not always can 
be obtained. Instead the numerical methods should be employed for this 
theory. Moreover, our discussion here will be restricted to the simplest case 
of decimation with the parameter rj = 2 on the tree level. Even such a simple 
transformation, however, can be of a use when the lattice calculations are 
concerned. 

Decimation diagrammatics once again consists of the internal propagator 

''(^^^E^iiF- '''' 

external legs 

Aik) = 2j2cos{ak^), (75) 



32 



inverse free decimated propagator G^^{K) and the vertices (Fig. 8, Fig. 9). 

The diagrams similar to one- dimensional (see previous Section), or equiv- 
alently matching conditions, determine form of the coefficient function for the 
effective action's quadratic part H^'^^K) ior rj — 2 as 

where 



8 (2 + cos(Xi/2) - cos(X2/2)) 8 (2 - cos(Xi/2) + cos(X2/2)) 

1 1 
^8{2 + Y:,cos{K,/2)) + 16 EMsm'(^M/4)j ^^^^ 

is an inverse decimated propagator. Up to the order 0{K'^) it can be ap- 
proximated as 



H(^)(X) -;^4sin^(i^,/2)-l5:(4sin^(X,/2))^ + |(;^4sin^(i^,/2))^ 

(77) 

On the other hand, this term should have the form (see Appendix C): 
H(^)(X) = - ^4sin^(X,/2) + (-^ + 2 c,) ^(4sin^(X,/2))V 

2c,{J2^sm'{K,/2)f. (78) 

Comparing these two expression, we find: 

C5 = 5/32; C6 = l/96. (79) 

Next, quartic terms in an effective action's expansion are calculated sim- 
ilarly to one dimensional case. Again, there are three contributions to these 
terms: classical term 

-(^/8) E(^x)^ 

X 

and two different diagrams (Fig. 13, Fig. 14). An effective theory, however, 
allows now for three different four derivative terms (see Appendix C) contrary 
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to the previous example where such term was unique. An effective action in 
momentum space has the form: 

-A^SfM^-^ J ^Hg),(X,L,M)7ri,7ri7rt,7rV^_^ (80) 



with 



U%{K, L, M) = Yi%,i{K, L, M) + c, Hj,.,,(i^, L, M) + c, H« L, M) 

*-ijkl 



+ c,Il'{K,L,M). (81) 



Once again the first term, 

nl-^iiK, L, M) = (^/4!) + L)) + (permutations)) , 

comes from expansion of the quadratic part of the covariant action in terms 
of pions. Next terms are: 

nl,i{K,L,M) = -(5/4!) L-4KE4sm'(i^M/2))(E4sin2(M^/2)) 

+ (permutations)) , 
H«.,,(X, L, M) = -(5/4!) (d.MY. 16 sin2(i^M/2) ^^n\Mj2)) 

+ (permutations)) , 

and 

Hj,,(i^,L,M) = -(5/4!) UmY^^HK,) sin(L,) sin(M^) 

xsin((-ii' - L - M)^))) + (permutations)) . 

To fully determine this quartic terms, we need now to fix three coefficients 
( cy, cs and cg). This can be done calculating the coefficient functions H^^\ 
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H*^^) for three linearly independent momenta configurations. f\ This will give 
us a system of three linear equations similar to eq. (^) for the coefficients. 
It is convenient to chose the following configurations: 

K=-L = M=ip, p), 

K={0, p), L = ip, -p), M = (p, 0), 

and 

K = L = M={p, 0). 
Then the coefficients C7, cs, cg will obey the equations: 

32/3 + 96 C7 + 48 cg + 96 cg = 19/6, 

245/48 + 78 C7 + 34 cg + 16 cg = 185/48, (82) 
17/2 + 120 C7 + 120 Cg - 72 Cg = 113/8 

with the solution: 

C7 = Cg = 0, Cg = -5/64 (83) 

The vanishing of C7 and cg is rather surprising. We do not see any obvious 
reason for this. 

On the tree level, one can still employ matching method to calculate next 
coefficients in expansion of an effective action: 

G^ivX, vY, ^Z) = G%{X, Y, Z). (84) 

Four point function as calculated in the original theory is given by the same 
diagram as in c? = 1 case. An effective theory, however, allows now for three 
different four derivative terms and thus there will be four contribution to the 
effective correlator Gf^li instead of two (Fig. 16). Thus to determine all the 
coefficients, we need to impose matching conditions for at least three inde- 
pendent momenta configurations. Numerical calculations gave the following 
result (a momenta configurations have been chosen the same as before): 

^Notice that technically decimation is nothing but summation over functions with an 
arguments shifted on aim. As a consequence, decimated expressions including functions 
like A(fc) can be sensitive to the sign of the argument even if the original expressions are 
not. Therefore, when calculating the decimated diagrams, one should take into account 
all different momenta permutations. 
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For P = 



-Q = R = (p, p); 



[|G<'')(P,Q,ii)|] = -^-l^ + ..., (85) 

GW(P,Q.fl) = -^ + ^iI±i£l±iii + ... (86) 
For P = (0,p), Q = (p, -p), ii = (p, 0); 

[|GW(P,Q,fl)|l = -i-H^ + ..., (87) 

GW(P. Q, R) = A+ ^ V '-'"'^ + ^ ^ + ... (88) 

For P = g = i? = (p, 0); 

[|G(^)(P,g,P)|] = -^-^ + ..., (89) 

G(4) (p, P) ^ -1^ + -^^/^ + (^Q/^) + (^Q/^) - « + ... (90) 

Matching conditions therefore give three hnear algebraic equations on the 
coefficients cy, Cg Cg and once again lead to the solution: 

C7 = Cg = 0, Cg = -5/64. (91) 



7 Conclusion 

To summarize we found a systematic way to perform RG of the decima- 
tion type m d > 1 perturbatively. We have seen during our discussion that 
perturbative decimation RG has rather complicated structure including con- 
siderable number of extra contributions compared to more customary ap- 
proaches (as momentum space RG) as well as some cumbersome numerical 
calculations. Here we would hke to some of its uses. 

First of all, the formalism we propose here, being based on the decima- 
tion RG transformations, pssesses all the advantages of this type of RG. As 
we have seen, it operates with original fields only and does not require any 
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(linear or nonlinear) transformations of variables. Besides, it preserves all 
the local relations including constraints. This in turn means that effective 
(coarse grained) theory will obey exactly the same local constraints as did 
original, and that no non-covariant terms will appear in the effective action. 
Among other applications, this opens a possibility to employ this formalism 
to study the critical phenomena, when strict control over symmetry proper- 
ties of model becomes particularly important. 

Another, compared to others P, ||], useful feature of proposed formalism 
is its perturbative character. This can provide us with systematic method of 
calculations in asymptotically free models and, what is even more essential, 
with a way to do controllable approximations. Hopefully, this side of proposed 
formalism will make it applicable in situations when such a control is essen- 
tial, as in above mentioned critical phenomena or in the recently proposed 
double strong- weak expansion approach It turns out that in asymptot- 



ically free theories there exist region in parameter space where both strong 
and weak coupling expansions are valid at the same time. Namely both the 
practical weak coupling a{g) = const g and strong coupling P{g) = const/ g 
expansion parameters are reasonably small. The "loop factors" l/(47r)^ in the 
practical weak coupling expansion parameter a{g) are partly responsible for 
this. In this scheme high frequency modes are integrated out perturbatively 
and the resulting effective action treated using strong coupling expansion. The 
symmetry preserving and controllable perturbative decimation technique is 
the most suitable tool for the first part of such calculations. Of course, to 
apply the method described here to the one dimensional models one should 
take into account correctly the nonperturbative instanton-like configurations, 
because in d = 1 the perturbation theory can be ill defined (one such exam- 
ple we considered in Section 5). In higher dimensions, however, the relative 
contribution of the nonperturbative configurations becomes less significant. 

We would like to stress also that the method described here unlike most 
of the other decimation (and exact RG in general) techniques, enables us 
to perform decimations not restricted to the simplest case rj = 2 only. The 
7] > 2 calculation just takes a bit more computer time. 
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8 Appendix A. d = 1 free fermion decimated 
action 



In this Appendix free Wilson fermions will be considered and decimated 
fermionic action in one dimension with rj — 2 will be derived via matching. 
The lattice action of the d — 1 Wilson fermions has the form: 



S = J2i'x^xy^y, (92) 



xy 



where 



^xy = ^[S{x -y + l)-5{x-y-l)] + am5{x - y) 



- r[5{x - y + 1) + 5(x - y + 1) - 25{x - y)\. (93) 
The Fourier transformed kernel is: 

QjJv 

A.{k) — ism{ka) + am + 2rsin^— . (94) 

Decimated propagator then reads: 

-1 ^ 1 

^ '^^^ " ^sin^^ + f + 2rsin2i^- ^^^^ 

First of all, notice that at r = this expression does not have the correct 
form of the fermion propagator, and is instead 

^-\^) = ,,Ti^. 2K - (96) 
m"^/4 + sm -j 

In particular, the corresponding effective action has no massless limit. The 
reason of such a strange behaviour is that actually in massless fermion theory 
without the Wilson term there is no way to build an exponential generating 
functional for decimated theory and therefore there is no way to define an 
effective action. Moreover, one can check that in such a theory correlators 
between the even (or odd) cites vanish, so that rj — 2 decimation in this case 
leads to the complete loss of information. 

In general case, the propagator takes the form: 

m^/4 + rm + (r"^ + l)sm -j + irsmK 
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that is the kernel of an effective action is 

MK) = — irsin K + -^—^ '- + (r + - sin^ • 98 

Notice that for massless theory effective action after field rescaling has 
the same form with the new Wilson parameter 

This means that the massless Wilson action in d = 1 is a perfect action and 
has a fixed point r = 1+. 



9 Appendix B. Diagrammatic derivation of 
d=l free boson decimated action 

Here we would like to reproduce the effective action in c? = 1 free boson model 
by means of the perturbation theory. Diagrammatics in this case consists of 
the internal propagator G{k), external leg A(fc) and the inverse decimated 
propagator A{K): 

^(^) = i^' A(fc) = ^cos(afc), A{K)^Asm'^. 
Diagrams contributing to the effective action are shown on Fig. 17 or 

-- o o +4 cn^-^- 



4 crv>~^ e e ■^-^-T^ 



Fig. 17. Diagrammatic representation of c? = 1 decimated action. 
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with 

^ 2 4 ..a^cos^aX,, 4 a^cosai^ ,,n , . r,K 

where denotes decimated expression: 

^TT ^ J-n/a ^^-^ 

The decimated blocs in an effective action are: 



ri a^cos{aK) \ ] 1 a 

CrW-^ O - LI 4sin22« IJ - 2' 



ri a^cos^{aK) n 1 a 

CrV>-^ ^-W^ - LI 4sin2 2|f IJ - 4sin2f 2" 



Thus, for the kernel of an effective action we obtain an expression 

A/r^N 2 4/1 a\ 4 / 1 a\\ . ^ . 

A(K) = h ^ Tw n \ Tw 4sm — = 4sm — , 

^ ^ a a2 \^4sin2f 2) o? 2) 2 2 ' 

(100) 

which clearly coincides with that obtained by the matching method. 



10 Appendix C. Quart ic part of the d = 2 a- 
model effective action 

In this Appendix an expression for the quartic part of an effective action of 
d = 2 (T-model will be given. This expression is necessary when the matching 
approach is employed. 

Decimation technique is applied to the unconstrained variables, "pions" 
TT* and therefore mainly gives us non-covariant quantities. On the other hand, 
an effective action is expressed in terms of constrained, covariant variables 
S"'. To reconstruct it from the non-covariant RG results, we need to re- 
express this effective action in terms of the fields tt* and then its coefficients 
can be identified by simple comparison. 
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The most general covariant effective action with up to four derivatives is 
given by an expression [21 

X 



1 



+C,{S''US''f + C8 Y.iS^'d.dp^f + C9 ^ 



(101) 

Here 5" (a = 1,2,..,A^) are 0{N) vectors on lattice normalized to unity: 
S*^ = 1, and we follow notations of with lattice spacing A = 1. 
To reconstruct coefficients C5, ..,Cg, one can solve the constraint: 



S^ = i^n\^l-g{7r^)^),t = l,..,N-l, 

and expand this action in terms of " pions" . To the fourth order in tt* and up 
to fourth derivatives, actions is: 



with quadratic part 

X . 

and quartic part 
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(102) 
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11 Appendix D. One loop contributions to 
the four derivatives terms of the effective 
action 

Here the leading in 1/N expansion one loop diagrams contributing to the 
four derivatives terms of the effective action are shown. All the diagrams can 
be divided into four groups: diagrams with all the external points coinciding 
(Fig. 18(a), 18(b)); with two pairs of the coinciding points(Fig. 18(c)); 
with two coinciding and two different points (Fig. 19) and with all different 
external points (Fig. 20). Diagrams 



Fig. 18. Contributions with 4 [figure (a)] and 2+2 (b) coinciding points. 

from different groups give different analytical expressions even ioi d = 1, 
where diagrams inside each group are proportional to each other. Thus, 
these groups can be calculated independently. This can provide us with 




o o 

o o 



(a) 



iP) 
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(a) (b) 
Fig. 19. Contributions of the 2+1+1 type. 



an additional consistency checks (cancellations of the negative powers of the 
external momenta inside the groups). 
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